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FOREWORD 


This report presents the analytical results of a study to define the 
diffuser requirements for sea level testing of high expansion ratio space 
engines. The work was performed by personnel In the Systems Analysis & Simu- 
lation Section of the Lockheed-Huntsville Research & Engineering Centrr under 
NASA contract NAS8-33981. The NASA Contracting Officer's Representative for 
this study was Mr. K, E. Riggs, EP23, 
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1, INTRODUCTION 


The design of a supersonic diffuser which will be compatible with a 
family of space engine systems and sea level test condl.i:lons Is a complex 
task. 


Lockheed-Huntsvllle has previously performed preliminary analyses of the 
diffuser requirements for sea level testing of the space engines, l.e«, the 
staged combustion and expander cycle engines. In the previous studies (Refs. 

1 and 2) specific diffuser designs were evolved. Since It was a design goal 
of MSEC to construct a "building block" space engine diffuser In which a whole 
series of space engines could be tested (for use on the Orblter Transfer 
Vehicle) this study effort was initiated. 

The study effort had three objectives: 

1. Review and update the analytical techniques used In diffuser 
design. 

2. Select the best of the analytical methods and conduct a 
parametric analysis of diffuser requirements. 

3. Apply the analytical techniques to a specific diffuser 
requirement. 

Diffuser starting transients and associated phenomena become more cri- 
tical as the nozzle area ratio Is Increased or the chamber pressure reduced. 

In fact, under the above conditions the diffuser start transients can either 
make or break a diffuser design. It Is fundamentally important to a success- 
ful diffuser design to thoroughly understand a particular diffuser starting 
phenomenon. The results of the analysis will also lead to operational seq- 
uences which must be performed. The "state of the art" of diffuser design 
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will be steadily pushed forward by the proposed MSFC Diffuser Facility and 
in many cases, the data will not exist in the literature. Therefore, extra- 
ploations of existing data using sound and fundamental engineering theory 
as a basis will be necessary. Another Important consideration for performing 
this analysis is the effect the physical presence of the diffuser will exert 
on the engine during the start transients. Ideally, the diffuser should simu- 
late the space vacuum environment. The physical presence of a diffuser, how- 
ever, will alter the engine starting', phenomena due to shock reflections, and 
start transient measurements may not be representative of a space vacuum start 
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2. TECHNICAL DISCUSSION 


2.1 LITERATURE REVIEW 

Research elfort has been directed In the past (Refs* 1 through 25 ) toward 
development of analytical methods for determining the performance of nozzle- 
diffuser systems. Analytical methods are avialable but are generally limited 
to the engine nozzle and entrance portion of the diffuser up to the second 
throat. In addition, empirical procedures have been developed which appear 
to be valid for predicting base pressure, overall diffuser performance and 
starting pressure ratio. The-iriany factors not Incorporated into the empir- 
ical methods prevent a purely theoretical approach to the diffuser design. 
Hence, past diffuser test experience is an Important factor in designing a 
new diffuser system. Some of the factors not considered in the empirical 
methods are: 

1. Heat transfer effects 

2. Initial boundary layer effects 

3. Real gas variable gamma effects 

4. Non-isoenergetic mixing effects 

5. Nonuniforra flow fields downstream of the nozzle 
throat includlug rf.gions of supersonic flow with 
bifurcated shocks and regions of subsonic flow. 

Hence the present "state of the art" for design and analysis of nozzle- 
diffuser systems consists of two basic procedures: 
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1* Empirical Methods ; These methods Involve calculations 
using available compressible flow equations using Isen- 
troplc, adiabatic and normal shock relations corrected 
by empirical factors which have been obtained experi- 
mentally for a variety of configurations* 

2. Theoretical Methods ; These methods utilize the computer 
to analyze the nozzle and diffuser flowfleld and are 
capable of treating one or more of the factors listed 
above which the empirical methods do not consider. 


2.2 DEVELOPMENT OP DIFFUSER DESIGN EMPIRICAL METHODS 

TH application of the empirical method used to design the engine diffuser 
is presented next. 


The general theory based on the simultaneous solution of the continuity* 
impulse, energy and state equations is presented in Ref. 3. The space engine 
parameters are presented in Table 3-1. As will be mentioned In section 3, 
the effective ratio of specific heats, y, for these engines is 1.22. This 
value will be used throughout rhls analysis. In order to reduce the diffuser 
entrance Mach number the diffuser entrance diameter is not made any larger 
than is absoutely necessary. Therefore the diffuser entrance diameter is 
generally fixed at 0.5 in. larger than the nozzle exit cutside diameter. The 
simultaneous solution of the continuity, impulse, energy and state equations 
results in an expression which may be expressed graphically as a function of 
y with Mach number as a parameter. This expression is: 







-Jl + ^M^ 


1 +yM^ 


where 

m = propellant mass flow rate in Ibm/sec 
2F = reaction force in Ibf 
R = gas constant in ft-lbf/lbm R 

T = total temperature in R 

** 2 
g^ ® gravitational constant in Ibm-ft/lbf-sec 
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M ■ Mach n jmber 
Y " ratio of specific heats 

A graph of this paraineter is presented In Figt 2-1* 

Straight Diffuser Analysis 

A sketch of this diffuser is shown in Fig* 2-2. The step-hy-step analysis 
follows. 

The diffuser entrance area to engine throat area ratio is 

Using this area ratio, the diffuser entrance Mach No. is obtained from 

isentropic tables aty® 1*22 

then from isentropic tables look up 

Pe/Pt* normal shock Pt 2 /Ptj^ and normal shock 

calculate the diffuser entrance force using the general thrust equation: 

F » Pj^ Ajjd + Y M^). 

Based on equilibrium chemistry data, the total temperature at the nozzle 
exit is obtained. Now compute the expression: 



Using this parameter in conjunction with Fig. 2-1 and assuming the dif- 
fuser L/D is approximately 8 or more, the exit Mach number is calculated. 
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Fig. 2-1 - Graphical Solution for Flow Parameter versus 
Mach Number, y « 1.22 
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Neglecting veil fcictlenf and beat transfer » ih«e! general thrust equation can 
be used to predict the diffuser exit pressure* R : 

0A 


P 


ex 


F 

*D <1 + V M./) 


As a check on this exit pressure we can use the normal shock pressure rise 
and compute P as 


ex 


P2/P1 


X P, 




X P 


GH* 


The value of P calculated by this method agrees favorably with the pre- 
vioMsly calculated value. The diffuse! starting pressure would be calcu-r 
lated as 


^CH start 


1,1 X P^jjj (14,7/Pgj^). 


Second Throat Diffuser Analysis 

The area contraction ratio is based on normal shock theory. The diffuser 
contraction ratio is computed thuslyj 

^ST 

•* Ag^ “ ^ ^ ^ST^^ 
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In this case, the exit force (see Fig, 2-3) is determined as 





in the absence of friction and heat transfer effects. Before we can proceed 
with the analysis further, a specific diffuser design must be established so 
that the forces acting on the two ramps can be calculated. In this study 
ramp angles were varied between 5 and 23 deg and flow properties were com- 
puted, The second ramp was added to increase the number of oblique shocks 
and therefore Increase the diffuser efficiency (Refs. 4 and 8), For this 
design the resulting exit force is computed as: 



\ 


- ®S 



) 


The flow parameter is then calculated, i.e,. 



From Fig, 2-1, the exit Mach number is read and the exit pressure is computed 
as 


ex 


ex 


*ex 


Now assuming a 60% subsonic diffuser efficiency at M the exit pressure, 

0X 

can be computed. 


P 


ex 
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The diffuser starting pressure Is then calculated as before, l.e.. 


^CH Start 


1.1 X Pjjjj X 14.7/Pgj^. 


The diffuser second throat contraction ratio can also be obtained from 
experimental data as shown In Fig. 2-4. As shown in Fig. 2-4, experimental 
data indicate that for a properly designed second throat diffuser, the con- 
traction ratio can be greater than one-dimensional normal shock theory suggests. 
The experimental curve will be used in the parametric analysis. 

Ejector Analysis 

An ejector is required when the diffuser exit total pressure is below 
16.17 psla (assuming 10% margin). The ejector must be designed to meet the 
requirements of ignition, full power and steady state performance at 10% power 
level. Ejector performance is enhanced when the ejector exit static pressure 
matches that of the engine flow at the diffuser second throat exit (Refs. 16, 

20, 24 and 25). In addition, ejector performance is further enhanced when 
the ejector flow expands radially outward (Ref. 16). The above statements 
infer that with fixed ejector flow rates, the ejector area ratio must be de- 
signed to be variable in order to match exit pressures at 10% and full power 
engine operation. In order to produce low cell pressures for ignition, a two- 
stage ejector is required with the first stage operating at a low flow rate 
(Fig, 3-4) and moderate Mach number. 

The ejector performance analysis method as outlined in Ref, i5 was used 
in all the parameter analysis work which will be discussed in Section 3. The 
second stage ejector was required to allow the first stage ejector to flow 
when the engine flow was not present, i.e., during start up and shut down. 

The philosophy used in the ejector analysis was based on the following: 

1. Analyses were conducted at the ejector exit and fully mixed 
stations without regard as to what transpired between (Ref. 

24). 
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Pig. 2-4 - Effect of Diffuser Inlet Mach Number on Maximum Contraction Ratio 
for Starting. (Ref. 13) 
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2. Primary and secondary flow are fixed to have the 
same static pressure at the inlet to the mixing 
section. 

3, Constant pressure mixing is assumed with subsequent 
flow area reduction down to the fully mixed region. 

4* The flow conditions before and after mixing are related 
to each other by the conservation requirements for 
mass, momentum and energy. 

5. Ejector length dimensions are based on available data 
to Insure fully mixed conditions. 

6. One-dimensional flow conditions with empirical corrections 
are assumed. 

The ejector analysis method outlined above allows one to relate the 
ejector performance specified by ejector pressure ratio, and ejector to 
engine mass flow ratio to the mixing section geometry given by its princi- 
pal cross-sectional areas. 

2.3 DEVELOPMENT OF DIFFUSER DESIGN THEORETICAL METHODS 

The General Interpolants Method (GIM) code (Ref . 2’7 ) was used to per- 
form the Navler-Stokes analysis of the diffuser flow field. A brief summary 
of the code follows. 

The GIM code employs a new methodology for constructing numerical 
analogs of the partial differential equations of continuum mechanics. A 
general formulation is provided which permits classical finite element 
methods and many of the finite difference methods to be derived directly. 

The GIM approach is new in the sense that it combines the best features 
of finite element and finite difference methods. The technique allows com- 
plex geometries to be handled in the finite element manner and operates on 
the integral form of the conservation laws. Solutions can be generated im- 
plicitly with the finite element analogs or by explicit finite difference 
analogs, which do not require a reduction of large systems of linear alge- 
braic equations (no matrix inverse) . A quasi-variatlonal procedure is used 
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to Introduce boundary conditions into the method and to provide a natural 
assembly sequence for combining the element equations Into the full domain 
equations. Attempts have been made In the literature to relate finite 
difference and finite element methods but have achieved limited success, 
and apparently no one other than Lockheed-Huntsville has used a combination 
of the two approaches In this unique fashion. 

As is the case with all attempts to solve partial differential equations 
by numerical approximations, the domain of Interest Is first discretized by 
appropriate subdivision into an assemblage of interconnected finite elements. 

A mesh generation is used in the GIM approach -hich incorporates general cur“ 
vlllnear coordinates, stretching transformations and bivariate blending to 
produce an automated mesh/element generation. Shape functions, based on a set 
of generalized Interpolants , are then chosen to describe the behavior over 
each element, We then proceed, as In the Method of Weighted Residuals by 
multiplying the descretlzed equations by a set of weight functions and inte- 
grating over the volume of the element. A quasl-varlatlonal procedure Is 
then used to construct the assembled system of equ,“:tlons from the element 
equations, and to Introduce boundary conditions Into the method. 

By choosiikg the weight functions equal to the shape functions, we reproduce 
the classical finite element nodal analogs. It Is at this point that we Intro- 
duce one of the important concepts of GIM; orthogonal weight/shape functions . 

By appropriately choosing the weight functions to be orthogonal to the shape 
functions, we can obtain explicit nodal analogs. Further, by a choice of 
arbitrary constants in the orthogonal weight functions, we can reproduce known 
finite difference nodal analogs, such as centered difference, upwind /downwind 
differences and the two-step MacCormak algorithm. As a result of this spatial 
discretization, we have reduced the partial differential equations to ordinary 
differential equations with '’time" as the independent variable. Any forward 
marching algorithm such as Euler, Runge-Kutta or predictor-corrector can be 
used to advance the sclution profiles in time. 
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The GIM formulation Is not a Finite Element method in the classical 
sense. Rather finite difference methods are used exclusively but the equa- 
tions are written in general orthogonal curvilinear coordinates. Transfor- 
mations are used to transform the physical planes into regions of unit cubes. 
The mesh is generated on this unit cube and the local metric coefficients 
generated. Each region of Che flow domain is likewise transformed and then 
blended to form the full flow domain. In order to treat completely arbi- 
trary geometric domains, different transformations may be employed in differ- 
ent regions. For this reason, we then transform the blended domain back to 
physical space. This allows the same set of equations to be solved in each 
region, with the local Jacobian of the transformation being the coefficients. 

This is not the classical treatment of many finite difference codes. To 
compute on the unit cube itself, one single transformation would be required 
for the entire domain. This would restrict the type of geometries that could 
be modeled. 

The following discussion summariz,es each module of the code. The GIM 
code is divided into four modules; (1) mesh generation/geometry; (2) nodal 
analog/matrix assembly; (3) unsteady integration; and (4) graphics. The mesh 
generation module accepts boundary ‘ geometry data, curve or line formula flags, 
and number of cuts in each coordinate direction. A set of general curvilinear 
maps is then used to subdivide each region into finite elements. Each region 
which is input is processed and then blended together. The output is a set 
of coordinates for each element along with the element coefficient matrices. 
The nodal analog assembly module takes the mesh data from a stored external 
file and performs via quasi-variational procedure, the assembly of the element 
equations into the full domain equations. At this point, the dynamic storage 
allocation is set up so that the unsteady integration module can integrate 
with virtually unlimited problem size. 

The unsteady integration module performs the actual computation of the 
flow by employing the boundary conditions selected by the user. The nodal 
analog at this point is arbitrary and any one of a number of schemes can be 
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selected depending on the problem being analyzed. The solution le marched 
forward In time for a specified number of steps or until a steady state is 
reached. The data display module reads the solution profiles from external 
storage (drum, tape) and prints, plots and maps the flow parameters. Fig. 

2~5 below is a block diagram illustrating the modular construction of the GIM 
code. 


Module I Modulo 2 Module 3 Module 4 



Fig. 2-5 - GIM Code Block Diagram 

The current version of the code can compute axisymmetric. and two- and 
three-dimensional viscous flows of an ideal gas in arbitrary geometric domains. 
The unsteady integration module is coded such that additional capability can 
be readily adapted such as different equation sets, other boundary values, 
virtually unlimited nodal points and time marching schemes. 

Through the choice of the arbitrary constants iti the orthogonal weight 
functions a forward-marching, or ’’quasi-parabolic” (QP) algorithm can be em- 
ployed in the GIM code. The forward marching technique requires only two com- 
putational planes to be in the computer at one time Instead of the entire flow 
field. This reduced storage greatly reduces computer runtime and is therefore 
advantageous for supersonic flows. 
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The diffuser flow field was modeled for the GIM code on the CGC CYBER 
203 computer at Langley Research Center. The entire flow field was modeled 
using the forward marching or QP algorithm. The QP solution resulted in a 
completely supersonic flow field with no normal shock in the tubular section. 
The results of the QP analysis at the entrance to the tubular section was to 
be input into a GIM code elliptic solution to compute the normal shock result- 
ing from Increased exit pressures. This elliptic solution was not performed 
due to program effort redirection by NA8A/MSPC. 
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3. PARAMETRIC ANALYSIS OP DIFFUSER REQUIREMENTS 


In this section of the study, a ^atametrie analysis was conducted using 
nozzle area ratio, e, chamber pressure, P , and degree of e.iector flow aug- 

a 

mentation, m^^^, as independent parameters. The diffuser configuration para- 
meters, such as area contraction ratio, Ag^/Aj^, length to diameter ratio, L/D, 
and the ejector design Mach number were design variables i At this time in the 
study, based on previous results (Refs. 1 and 2), it was decided to fix the 
ejector inlet total conditions to correspond closely to the conditions which 
exist at the turbine exit in a turbojet engine. The ejector inlets were 
thus fixed, but the mass flow to each ejector was varied as required. To 
accomplish this study the one-dimensional empirical methods presented in 
Section 2 of this report were programmed on the in-house PDP-11 digital com- 
puter. Then, using the results generated, diffuser commonality and "building 
block" approach was considered. In the above analysis, a 10% starting margin 
(based on chamber pressure) was used. 

Table 3-1 presents the characterltistlcs of five candidate space engines. 
Since the "building block" diffuser must accommodate all of these engines, the 
diffuser duct inlet diameter D^^ was set at 73.5 in. 

The thermochemical properties for the LOX-LH 2 propellant system were 
obtained from Ref. 26 and previous chemical equilibrium combustion computer 
program runs made for similar operating conditions. The results of this data, 
indicate that the ratio of specific heats, y» varies from 1.14 in the com- 
bustion chamber to 1.32 near the nozzle exit plane. It was determined that 
for use in the simplified one-dimensional flow diffuser analysis that an 
effective y of 1.22 would yeild satisfactory nozzle exit flow conditions. 
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Table 3-1 - FIVE CANDIDATE SPACE ENGINES FOR THE OTV 


Parameter 

Thrust* Full 
Thrust, Low 

Maximum Test Duration @ MR « 6.0 
Full Thrust 
Low Thrust 
Glmbal Capability 
Propellants 

Mixture Ratio, Full Thrust 
Low Thrust 
Nozzle Area Ratio 
Engine Evelope: 

Outside Diameter @ Noz, Exit 
Inside Diameter A Noz. Exit 
Length, Glmbal Pad to Noz. Exit 
Length, Glmbal Pad to Inlet Flange 
LOX 
LH2 

Engine Weight 

Chamber Pressure, Full Thrust 
Cahmber Pressure, Low Thrust 
Nozv Exit Wall Press., Full Thrust 
Noz. Exit Wall Press., Low Thrust 
Total Flow Pate, Full Thrust 
Total Flow Rate, Low Thrust 


Units 

AECE-A 

AECE-P 

AECE-R 

ASE RLIO-IIB 

lb 

15000 

15000 

15000 

2000 

15000 

lb 

2000 

1500 

1800 

1850 

1500 

sec 

1200 

1200 

1200 

1200 

1200 

sec 

2500 

2500 

2500 

2500 

2500 

— 

none 

none 

none 

none 

none 


L 0 X/LH 2 

LOX/LH^ 

L 0 X/LH 2 

LOX/LH 2 

LOX/LR, 

— 

6.0 

6.0 

6.0 

6.0 

6.0 


6.0 

6.0 

4.0 

2,0 

6.0 

— 

473 

642 

625 

400 

205 

In. 

62.7 

66.1 

63.25 

58.08 

73.0 

In. 

60.7 

64,1 

61.25 

56.08 

71.0 

In. 

120 

114 

117 

100 

110 

In. 

12 

12 

12 

27.1 

10 

In* 

15 

15 

15 

36.87 

10 

lb 






Psla 

1200 

1505 

1539 

2028 

400 

Psla 

160 

150 

198 

187 

40 

Psla 

0.196 

0.163 

0.172 

0.406 

0.19 

Psla 

0.026 

0.016 

C.022 

0.037 

0.019 

Ib/sec 31.4 

31.2 

31.2 

42.01 

32.6 

Ib/sec 4.2 

3.2 

3.7 

4.06 

3.2 
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No Y changes wave considered In the analysis. Additional study effort 
should be undertaken to determine the effects of the real gas variable y 
on the diffuser design* The ratio of specific heats used for tha ejector 
flow was Y** 1.33. 

Basically two specific tasks were undertaken and the results are pre- 
sented In each one. 

3.1 PARAMETRIC ANALYSIS FOR AN UNAUGMENTED DIFFUSER 

In this task, both high chamber pressure (high thrust) and lew chamber 
pressure (low thrust) operation were analyzed as to the maximum allowable 
nozzle area ratio which tested In an unagumented diffuser. The results are 
presented In Figs. 3-1 and 3-2, Figure 3-1 presents the results for the low chamber 

pressure operation, and Fig. 3-2 presents the results for the higher chamber 
pressure operation. As shown In Fig. 3-2 only the ASE engine can be run in 
a properly designed unaugmented diffuser. 

3.2 PARAMETRIC ANALYSIS FOR AN AUGMENTED DIFFUSER 

For the performance ot this task the original computer program developed 
to accomplish Task 3*1 was modified to analyze diffuser and ejector perfor- 
mance together and ejector performance alone with no engine flowing so that 
ejector blank off performance could be evaluated. The computer program has 
the capability to add up to two ejectors (working in stages) as required to 
make the engine diffuser combination flow against ambient pressure. A +10% 
margin on chamber pressure Is Incorporated Into the diffuser design computer 
code (DIFEJJ) . 

For this task a decision had to be made on a choice of ejector 
fluid. As mentioned earlier in this section, turbojet exhaust was selected 
over steam. The reasons for this selection were: 
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Chamber Pressure, P (psia) 


I’ig. 3-1 - Unagumented Diffuser Limiting Area Ratios for Various Space 
Engine Chamber Pressures (Low Thrust) 
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Fig, 3-2 - Unaugmented Diffuser Limiting Area Ratios for Various Space 
Engine Chamber Pressures (High Thrust) 
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1. Engine test durations o£ 30 min are required, consuming 
ejector fluid at a rate in excess of 300 Ib/sec. An equiva- 
lent steam ejector system would require a building 100 ft 

by 300 ft by three stories high just for the steam generating 
equipment. 

2. Quick startup capability using turbojets 

3. More reliable operation than a steam system, and 

4. Lower capital Investment and much lower operating cost. 

Turbojet engine literature was reviewed as to the exhaust flow properties. 
Typical turbine exist flow properties are presented in Fig, 3-3 which shows the 
"building block" diffuser and ejector configuration. One main requirement 
placed on the ejector system was that it be able to evacuate the diffuser 
down to at least 0,4 psia. Fig, 3-4 shows the typical variation of diffuser 
evacuation pressure with ejector flow rate for two different ejector design 
Mach numbers. As a result of studies like these It was decided to fix the 
ejector flow rates at 32,6 Ibm/sec for ejector 1 and 293.4 Ibm/sec for ejector 
2. These flow rates correspond to several available turbojet engines such as 
the General Electric J85 and the Pratt & Whitney JT-llB-2 (J58) . 

The results of this parametric study are shown in Fig. 3-5, The results 
indicate that with a properly designed diffuser and two stage ejector system, 
space engines with area ratios up to 2000:1 and chamber pressures down to 40 
psia can be tested. For all the envelope shown, the diffuser evacuation pres- 
sure was between 2,0 mmHg and 4,0 mmllg. The first ejector Mach number varied 
from 1,2 to 4.54, The second ejector Mach number varied from 2,44 to 2.66. 

The operation of the ejector/dif fuser system is as follows: 

1. The specific space engine is located on the operational 
envelope (Fig. 3-5) in terms of engine chamber pressure 
and area ratio. The space engine will fall in either one 
of three categories. 

• No augmentation required for space engine operation 

• First stage ejector required for space engine operation 

• Two-stage ejector required for operation. 
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2« If the space engine requires only the first-stage ejector 
for operation* both stages must be used during start up 
and shut down* The system also features a gaseous nitro- 
gen (GN 2 ) purge system as an extra safety precaution to 
eliminate the possibility of detonation on start up due 
to fuel lead. Prior to startup the two-stage ejector 
is designed to evacuate the diffuser dowi\ to 0.06 psia 
or 3 mniHg, An equivalent steam ejector system would re- 
quire a three stage ejector system. 
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4. DIFFUSER REQUIREMENTS EXAMPLE 


In this study, the general rocket exhaust diffuser requirements were in- 
vestigated for the sea level testing of the RLIO-IIB (Table 3-1) Pratt Whitney 
Expansion Cycle engine. The operating parameters for the engine are presented 
in Table 3-1. A preliminary engine exhaust diffuser design ws developed and 
is presented in this section. Detailed engine nozzle and diffuser flowfleld 
analyses were performed. The RUO-llB engine diffuser design was developed 
first since this engine will be available first and its diffuser operating 
data will be invaluable in establishing the diffuser requirements for the 
other space engines. The results obtained from this preliminary study are 
summarized in this section. 

4.1 GEOMETRY DESIGN DETAILS 

The diffuser design function must be performed before detailed flow 
analysis can be performed. In order to maintain the required low nozzle exit 
pressures for the nozzle to flow full, and still exhaust to atmospheric con- 
ditiontSj the supersonic diffuser must exhibit a pressure recovery better than 
normal shock recovery based on a cavity free stream Mach number of approxi- 
mately 6. Based on a review of the literature. Ref. 1 through 25 and various 
LMSG reports this requirement necessitates the use of a specially designed 
second throat. The second throat minimum area is still based on the normal 
shock method since during starting, the second throat area must still be large 
enough to pass the normal shock flow. A range of second throat inlet ramp 
angles between 5 and 23 deg was investigated in this study. To obtain good 
total pressure recovery in as short a distance as possible, a double ramp 
second throat was selected (Ref. 4), This type of design was built and tested 
in the Ref, 4 study with excellent results. One interesting flow phenomena 
situation which occurred during this analysis and was borne out in the Ref, 4 
tests, is the length of the cylindrical section between the first and second 
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ramps as shovm in Fig, 4-1. The analysis indicated that if this cylindrical 
length is too long, the flow will expand back to a high Mach number with sub- 
sequently larger total pressure losses. This phenomenon was demonstrated 
in the Ref, 4 study when the diffuser starting pressures were reduced substan- 
tially when the length-to-diamter ratio of the cylindrical section was reduced 
from 1,53 to 0.51, 

4,2 ONE-DIMENSIONAL ANALYSIS RESULTS 

The overall diffuser design (Figs. 4-1 and 4-2) is based on the experi- 
mental and analytical results presented in Refs. 3 and 4 and discussed further 
in Section 2. The overall diffuser design was calcualted using the empirical 
one-dimensional diffuser/ejector design code (DIFEJJ) developed in this study. 

The first stage ejector exit conditions shown in Figuco 4-3 for this design 

a r e a s £o Hows , 1 . e , , 


• High Thrust operation (Pc * 400 psia) 

Ejector Mach number « 1.72 
Blank-off pressure = 0.068 psia 
Flow Rate « 32.6 Ibm/sec 

• Low Thrust operation (Pc « 40 psia) 

Ejector Mach number =3,15 ) 

Blank-off pressure = 0.055 psia | Not shown in Fig. 4-3 
Flow rate = 32,6 Ibm/sec ) 


The second stage ejector operating conditions are shown in Fig. 4-4. The second 
stage ejector Mach number for the high thrust and low thrust operation are 2,52 
and 2,39, respectively. The high thrust operation variables are presented in 
Figs. 4-3 and 4-4. 

The high thrust mode of operation will impose the most severe heat loads, 
therefore the flowfield details were developed for it. 
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4.3 DETAILED ANALYSIS RESULTS 

In this section the results of the detailed analysis for the diffuser 
presented In Figs. 4»1 through 4-4 are presented. The results presented are 
for the flov^leld parameters, heat loads, water cooling requirements, afterf- 
burning in the ejector mixing section and startup transient concerns. 

4.3.1 Flow Field 

The diffuser flow field was computed using the GIM code, discussed In 
Section 2.3, in the axlsyraraetrlc forward-marching mode. The geometry used 
in the analysis is shown in Fig. 4-5. The geometry was Initiated downstream 
of the nozzle throat In order to ensure supersonic flow throughout the flow 
field. The nozzle was included In the analysis In order to obtain a good 
flow profile at the nozzle exit. The nozzle throat area used to compute the 
start line condltons was adjusted to accommodate the constant gamma assumption 
(Y “ 1*22) . Chamber conditions of 400 psia pressure and 5800 R temperature 
were used to compute the start line conditions along with a fluid flow rate 
of 32.6 Ibm/sec and a mixture ratio c.f 6.0. The nozzle configuration was 
obtained from Pratt & Whitney through NASA-MSFC. 

The diffuser geometry consisted of 364 planes containing 25 nodes each for 
a total of 9100 nodes. A free-sllp tangency boundary condition was Imposed on 
the diffuser wall. If the elliptic analysis of the tubular section had been 
performed, no-slip boundary conditions would have been imposed in order to 
develop the boundary layer in the tube. 

Pressure contours for the resultant flowfleld analysis are shown in Fig. 
4-6. The pressure contours are not normalized (P^ “ 1,0) and represent pres- 
sures in psfa. The wavy contours near the nozzle axis are due to plotter dis- 
cretization and depict regions of nearly constant pressure. A strong oblique 
shock is generated when the nozzle flow encounters the diffuser wall. A 
second shock occurs when the diffuser necks down at a 6. deg angle. This 
second shock blends with the first as they reach the axis, A third weak ob- 
lique shock is generated at the second 6 deg ramp, 
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The interaction of these shocks and their respected reflected shocks, 
along with the compression on the six~degree ramp, creates a very complex 
n.ow pattern just prior to entering the tubular section. The stronger re- 
flected shock continues into the tubular section with decreasing strength. 

The flow conditions along the diffuser wall, required for heat transfer 
analysis, are presented in Table 4-1. 

4.3.2 Diffuser Heat Load Analysis 

Heating rates to the inside wall of the diffuser were calculated using 
the Lockheed Multiple Pressure Gradient Program of Ref. 28. All heating rate 
calculations were based on the 100% thrust flow field assuming fully turbulent 
flow. The governing equations for the analysis are presented in Appendix A 
and were extracted from Ref. 28 for the convenience of the reader. 

Thermal and transport properties used by the program are calculated using 
a model based on the Lennard-Jones potential Intermolecular force model as 
documented in Ref. 28. The mass-fractions of the chemical species used in the 
analyses were determined using the NASA-Lewis CEC program. Chemical equil- 
ibrium was assumed at all locations. Chemical species considered in the 
heating analyses include H, H 2 , HO 2 , H 2 O, 0, OH, 02» Heating rates based on 
a cold wall (540 R) are shown on Fig. 4-7 as a function of axial distance along 
the length of the diffuser, A peak in the heating occurs at a point just down- 
stream of the point where the flow first impinges on the diffuser wall. The 
maximum heating value is approximately 87 Btu/ft -sec which occurs in the 
psuedo-shock flow region of the second throat tubular section. The change in 
slope of the heating versus distance curve at appoximately 22 in. is due to the 
flow expansion caused by the diffuser geometry changing from a conical to a 
cylindrical section. 

Figure 4-8 shows the gas temperature. Fig. 4-9 shows the gas pressure, 
and Fig. 4-10 shows the gas Mach number along the diffuser wall. 
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Table 4-1 - DIFFUSER WMiL CONDITIONS (Pc - 400) 


X (ft) 


^wall 


(ft) 


Density 

(Ibin/ft^) 


Static Press* Mach 

(psfa) No. 


e.46H 

3.0620 

0.000194 

39.760 

4.405 

0.583S 

3.0620 

0.000226 

48,890 

4.257 

0.6856 

3.0620 

0.000256 

57,660 

4. 146 

0.?B?? 

3.0620 

0.000282 

65.340 

4 064 

0.88f8 

3.0620 

0.000303 

71.760 

4.011 

0.99S8 

3.0620 

0.000318 

76.850 

3.953 

1.0949 

3.0620 

0,000352 

89.070 

3.807 

(First Ramp) 





1 . 1850 

3.0530 

0.000383 

99.180 

3.744 

1.2750 

3.0430 

0.000411 

108.480 

3.691 

1.3660 

3.0340 

0.000434 

116. 160 

3.650 

1.4560 

3.0240 

0.000451 

122.030 

3.618 

1.5470 

3.0150 

0.000462 

126.010 

3.595 

1.6380 

3.0050 

0.000466 

128. 150 

3.579 

1.7280 

2.9960 

0.000466 

128.530 

3.520 

1.8190 

2.9860 

0.000461 

127.330 

3.568 

1.9090 

2.9770 

0.000452 

124.760 

3.572 

2.0000 

2.9670 

0.000440 

121.090 

3.581 

2.0910 

2.9580 

0.000425 

116.640 

3.595 

2. 1810 

2.9480 

0.000409 

111.730 

3.612 

2.2720 

2.9390 

0.000393. 

106.640 

3.633 

2.3630 

2.9290 

0.000348 

95.080 

3.626 

(Straight Sec 

tion) 




2.4520 

2.9290 

0.000302 

80.960 

3.697 

2,5410 

2.9290 

0.000264 

69.490 

3.770 

2.6310 

2.9290 

0.000228 

58.790 

3.852 

2.7500 

2.9290 

0.000205 

51.870 

3.915 

2.8700 

2. 9290 

0.000190 

47.550 

3.959 

2.8980 

2.9290 

0.000181 

44.870 

3.988 

3. 1080 

2 . 9290 

0.000175 

43. 160 

4.007 

3.2280 

2.9290 

0.000170 

41.750 

4.023 

3.3770 

2.9290 

0.000166 

40.690 

4.035 

3.5260 

2.9290 

0.000 163 

39.820 

4,645 

3.6760 

2.9290 

0.000161 

39.020 

4.056 

3.8240 

2.9290 

0.000157 

38.270 

4.066 
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Table 4-1 - DIFFUSER WALL CONDITIONS (Pc - 400) Cont’d 


X (ft) 


Density 
(Ibm/f t^) 

Static Press, 
(psfa) 

Mach 

No. 

3.9?30 

2.9290 

B.00015S 

37.540 

4.076 

4. 1530 

2.9290 

0,000153 

36.830 

4,086 

4.2F20 

2.9290 

0.000150 

36. 130 

4.097 

4.4210 

2.9290 

0.000140 

35.440 

4.107 

4.5700 

2.9290 

0.000146 

34.770 

4.117 

4.F190 

2.9290 

0,000144 

34.240 

4. 125 

4. 8390 

2.9290 

0.000142 

33.720 

4. 133 

4.9580 

2,9290 

0.000140 

33.210 

4. 141 

5.0770 

2.9290 

0.000139 

32.830 

4. 147 

5.1G70 

2,9290 

0.000137 

32.460 

4. 153 

5.2560 

2.9290 

0.000136 

32. 100 

4. 159 

5.3460 

2.9290 

0.000149 

32.890 

4.008 

(Second Rami)) 





5.4500 

2.9180 

0.800163 

41.130 

3.946 

5.5550 

2.9070 

0.000174 

44.780 

3,896 

5.G590 

2.8960 

8,008183 

47.710 

3.859 

5.7640 

2.8850 

0.000190 

49.920 

3.832 

5.8630 

2.8740 

0.000194 

5 1 . 450 

3.813 

5.9740 

2,8630 

0.000197 

52.370 

3.801 

6.0800 

2.8520 

0.000198 

52.780 

3.795 

6. I860 

2. 84 IP 

0.000198 

52,760 

3.793 

G . ..-.’riO 

2.8300 

0.000197 

52.400 

3.794 

6.4000 

2.8180 

0.000194 

51.800 

3.798 

6 . 5070 

2.8070 

0.000192 

51.030 

3.805 

6. 6 160 

2,7960 

0.000109 

50.160 

3.812 

6.7250 

2.7840 

0.000186 

49.260 

3.821 

6 . 8340 

2,7730 

0.000183 

48,360 

3.829 

6.9450 

2.7610 

0.000180 

47.490 

3.838 

7,0570 

2.7490 

0.000176 

46 . 680 

3.846 

7 . 1690 

2.7380 

0.00017S 

45.930 

3.853 

7.2830 

2.7260 

0.000 IP'S 

45,240 

3,860 

7.3980 

2.7130 

0.000171 

44.620 

3.866 

7.5140 

2.7010 

0.000169 

44,050 

3.872 

7.6310 

2.6890 

0.000167 

43.530 

3.877 

7.7500 

2.6760 

0,008165 

43.050 

3.882 

7.8710 

2.6640 

0.000164 

42.600 

3.887 

7.9920 

2.6510 

0.000163 

42. 170 

3.891 

8 . 1 tnO 

2.6300 

0.000161 

41 .770 

3.895 

8.2420 

2.h2Sa 

0.000160 

41.390 

3.899 

8.3690 

2.6110 

0.000159 

41.020 

3.903 

8.4980 

2.5980 

0.000157 

40.670 

3.987 

8.6300 

2.5840 

0.000156 

40.330 

3.911 

0,7640 

2.5700 

0.000155 

40.000 

3.914 

8.9000 

2.5560 

0.000154 

39.700 

3.918 

9.0390 

2.5410 

0.000153 

39.410 

3.921 

9.1800 

2.5260 

0.000152 

39 . 130 

3.924 

9. 3250 

2.5110 

0.000151 

38.880 

3.926 
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Table A-1 - DIFFUSER WALL CONDITIONS (Pc - 400) Concl'd 


S (ft) 


^wall 


(ft) 


Density 

(Ibm/ftS) 


Static Press. Mach 

(psfa) No. 


9.4720 

2,4950 

0.000150 

38.640 

3.929 

9. S 230 

2.4800 

0.000130 

38.430 

3.931 

9,7770 

2.4630 

0.0 U 0149 

38.240 

3.933 

9.9310 

2.4470 

0.000148 

38.080 

3.934 

10. OB 00 

2.4310 

0.000148 

37.960 

3.935 

10.2300 

2.4160 

0.000148 

37.860 

3.936 

10.3700 

2.4010 

0.000147 

37.780 

3.936 

10.5100 

2.3860 

0.000147 

37.730 

3.936 

10.6500 

2.3710 

0.000147 

37.700 

3.936 

10.7900 

2.3570 

0.000147 

37.690 

3.935 

10.9200 

2.3430 

0.000147 

37.700 

3.934 

11.0600 

2,3290 

0.000147 

37.720 

3.934 

11.1000 

2.3150 

0.000147 

37.770 

3.932 

11.3100 

2.3020 

0.000147 

37.840 

3.931 

11.4400 

2.2890 

0.000147 

37.930 

3.929 

11.5600 

2.2760 

0.000148 

38.010 

3.927 

11.6800 

2.2630 

0.000140 

38. 130 

3.925 

11.0000 

2.2506 

0.000148 

38.260 

3.923 

11.9200 

2.2380 

0.000149 

38.400 

3.920 

12.0400 

2.2260 

0.000149 

38.560 

3.918 

12.1600 

2.2130 

0.000150 

38.740 

3.915 

12.2700 

2.2010 

0.000150 

38.930 

3.912 

12.3000 

2.1890 

0.000151 

39.130 

3.908 

12.5000 

2. 1780 

0.000152 

39.360 

3.905 

12.6100 

2. 1660 

0.000152 

39.600 

3.901 

12.7200 

2.1540 

6.000153 

39.858 

3.897 

12.8300 

2. 1430 

0.000154 

40. 130 

3.893 

12.9400 

2. 1310 

6.000155 

40.420 

3.886 

13.0600 

2.1200 

0.000156 

40.740 

3.884 

13.2600 

2.0970 

0.000158 

41.450 

3.874 

13.3700 

2.0860 

0.000159 

41.840 

3.868 

13.4700 

2.0750 

0.006161 

42.280 

3.862 

13.5800 

2.0640 

0.000162 

42.760 

3.855 

13.7900 

2.0420 

0.000166 

43.900 

3.840 

13.8900 

2.0310 

0.080168 

44.590 

3.831 

14.8000 

2.0200 

0.000170 

45.380 

3.821 

14. 1000 

2.0090 

0.000173 

46.320 

3.808 

14.2100 

1.9980 

0.000175 

47.260 

3.895 

26.2100 

1.9980 

0.000175 

47.260 

3.895 

38. 1900 

1 .9980 

0.001756 

1087.200 

0.351 


(End of Second Ramp) 
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Axial Distance - Inches 

Fig. 4-8 Diffuser Local Temperature Distribution 
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4.3.3 Dl££user Cooling Water Requirements 

The diffuser was divided into six separate segments to perform the 
cooling water analysis. The separate segments are shown in Fig. 4-7 . The 
results of this analysis are shown in Table 4-2. The maximum gas side fiJr: 
coefficient was used over the entire section. The water flow rates were de- 
signed to be the same through each section with the water bulk temperature 
rise being a variable as shown in the last column of Table 4-2. In this way, 
the same water can be used from one end to the other, starting at the ejector 
end. This design minimizes the water flow requirements with the total bulk 
temperature rise being only 71 F. The diffuser wall thickness between the 
hot gas side and the water was selected at 0.125 in. The wall would be struc- 
turally stiffened by the water channel ribs and the outer water channel skin. 
The entire diffuser could be made from mild; steel sheet since the loaximum gas 
wall temperature is only 442 F. As a result of this analysis, no thermal 
problems are anticipated with the diffuser. 

4.3.4 Afterburning Analysis 

The engine flow exiting from the diffuser tube into the first ejector con- 
tains excess H 2 gas. The ejector gas contains excess hot air. The purpose of 
this analysis was to determine if any appreciable afterburning would occur in 
the first ejector mixing tube thereby raising the temperature and possibly 
choking the flow. The pertinent flow variables prior to mixing are shown in 
Fig. 4-3. The central gas core passes a nomal shock wave prior to its exit 
from the tube, thereby raising its temperature up to 5700 R. Also shown for 
comparison purposes in Fig. 4-3, are the results that were obtained using the 
area averaged GIM code upstream flow results and processing the flow across a 
normal shock. Excellent agreement (within 8% on pressure) is obtained between 
the 1-D DIFEJJ CODE and the GIM code for the main engine flow exiting the 
diffuser tube. 

The laser And Mixing Program (LAMP) Ref. 29, was used to perform the 
afterburning analysis. Tiiis program has been used in the past by Lockheed to 
study the Space Shuttle main engine (SSME) exhaust/ambient air mixing phenomena 
and good comparisons were obtained with the available measured data. Since 
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the SSME uses the identical propellant system as the R1.10‘'>IIB, the same 
finite rate chemistry package was used in this analysis. 

The results obtained are subject to the following assumptions: 

1. Mixing and reactions proceed at constant pressure. 

2 , Due to the current setup of the XMP code for axisymmetric 
mixing, conditions on the outermost streamline (i.e,. the 
outer edge of the ejector flow) will remain constant. This 
implies a constant ejector wall temperature, as well as con- 
stant velocity and gas composition along the wall. 

The results indicated that nc appreciable afterburning will occur 
and that the effect of the ejector flow is tc c<-'>l the hot diffuser gas with 
a subsequent flow area reduction. 

4.3.5 Engine/Diffuser Startup Transients 

During the RLIO-IIB engine startup transient which takes about 3 sec. 
the H 2 fuel leads the oxidizer. The fuel flow builds up to a maximum of 5.0 
Ibm/sec in 3 sec duration. During this time the igniter is firing. Combus- 
tion will not occur until a minimum of 4 psla chamber pressure is reached in 
the presence of some oxidizer. The concern is what happens when the leading 
H 2 fuel reaches the relatively ^jarra turbojet oxygen-rich exhaust. This problem 
is shown pictorlally in Fig. 4-11. As can be seen in the figure, the ejector 
blank-off conditions in the diffuser tube are pressure « 0.068 psia and 
teroncrature of ^ 29R. This condition will not promote detonation of the 
H 2 fuel. The H 2 fuel will be mixed with the ejector flow and swept downstream 
and burned. To abjtolutely preclude any possibility of detonation, a purge flow 
of N 2 has been added as shown in Fig. 4-11 and previously in Fig. 3-3. 
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5. CONCLUSIONS 


As a result of this contract, a very useful diffuser /ejector design com- 
puter program was developed. Using empirically modified one-dimensional flow 
methods the diffuser ejector geometry is specified by the code. Results ob- 
tained using the detailed GIM code verified the design code results for cal- 
culations up to the end of the diffuser second throat. 

This study has demonstrated that through careful extrapolation of existing 
analytical methods and experimental data, a two-stage diffuser/ejector building 
block system can be designed which will accommodate any potential OTV space 
engine . 

The results of this study have demonstrated that through careful design, 
an ejector system using two commonly available turbojet engines feeding two 
variable area ratio ejectors is indeed feasible. In fact, the use of turbo- 
jet engines as ejector drivers is much more practical than the use of steam. 

This study has defined diffuser requirements for sea level testing of 
high expansion ratio space engines. The analysis has demonstrated numerically 
the practicality and feasibility of a "building block" engine diffuser/ejector 
system. The next step would be to define and evolve the design further and 
study the possibility of building a prototype system. However, due to scale 
model boundary layer effects, a prototype system smaller than 1/4 scale may 
not provide meaningful results* 

Future analytical effort should involve additional GIM code analysis of the 
two stage ejector mixing process. 
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NOMENCLATURE FOR APPENDIX A 

Symbols 



specific heat at constant pressure, Btu/lb*^R 

8h 

enthalpy conductance, Ibm/ft^-sec 


adiabatic wall enthalpy, Btu/lb 

J 

conversion factor, 778 ft-lb/Btu 

Nu 

Nusselt number, dimensionless 

Pr 

Prandtl number, dimensionless 

q 

convective heat transfer rate, Btu/sec-ft^ 

r 

recovery factor, dimensionless 

Re 

Reynolds number, dimensionless 

S 

distance along a flow streamline, ft 

T 

temperature, °R 

V 

velocity, ft/sec 


flow characteristic length, ft 

Greek 

e 

shape factor, ft 

n 

cone half angle, deg 

p 

density, Ib/ft^ 

Subscripts 

e 

property evaluated at the edge of the boundary 
layer 

r 

property evaluated at the defined reference 
condition 

w 

property evaluated at the wall 
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Appendix A 


In this Appendix, the boundary layer heat transfer theories used in the 
diffuser heat load analysis are presented. 

• Laminar Boundary Layer 

Heat transfer through a laminar boundary layer is calculated using the 
integral form of the energy equation. Effects of variable freestreara velocity, 
density, and pressure are accounted for through use of appropriate trans- 
forms of the flat plate solution. Non-constant properties through the boundary 
layer are also accounted for. 

In Ref. A-1 (among others) it is shown that the energy equation can be re- 
duced to an equation with enthalpy as the dependent variable under a number 
of different conditions Including some cases where there is dissociation, chemi- 
cal reaction and mass diffusion within the boundary layer. The complete dev- 
elopment of the enthalpy form of the equation and the evaluation of the sur- 
face heat rate are presented in Ref. A-1. 

Based on examination of several "exact" laminar boundary layer solutions, 
Eckert (Ref. A-2) recommends that the effects of variable gas properties 
through the boundary layer can be accounted for by simply evaluating the pro- 
perties at a "reference enthalpy" and using these values in the constant 
property solutions as obtained by Blasius. Based on thif, method, the convec- 
tive heat rate to the wall is evaluated vsing 
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where the adiabatic wall or recovery enthalpy is 


whe re 



,v* 

c 

2Jg 


t ?= 



(A-2» 


(A- 3) 


The enthalpy conductance is evaluated from the Stanton number (St) 


St* = 


8v 


p*v; 


Nu* 


Re* Pr* 


(A-4> 


The starred (*) properties refer to properties evaluated at a temperature 
corresponding to a ’'reference enthalpy," h* , where 


h* sh^ + 0.5 (h^ - hg) + 0.22 (h^ - h^) . (A- 5) 


It can be readily seen that an iterative solution is required for Eqs. (A-2), 
(A- 3) and (A-4). 

The Blasius flat plate constant-property boundary layer solution yields 



0.332 Pr 


1/3 


Substituting Eq. (A- 6) into (A-4) yields 



(A-6) 


(A-7) 
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Substituting, Rc « 


P v^Xl 


and Eq. (3.80) into Eq. yields 



0.332 

7^^ 


X 

*L 





(A-8) 


The characteristic length (X|^) as used in Eq. (A-8) is obtained for variable 
property flow by numerically integrating the following equation along a flow- 
field streamline, (Refs. A-1 through A- 3). 


X 


L 


* * ,, 
p ui y. 


f p* 


e ^ dS. 


(A** 9) 


The parameter S is the wetted length along the streamline. For a flat plate 
the parameter e is set equal to a constant and thereby divided out, and for 
constant inviscid flow properties, the characteristic length (X^) be ‘.omes 
equal to the flat plate length. For axisymmetric flow about a cone with con- 
stant inviscid properties, the equivalent length, X^^, becomes 



For a cone with a cone half angle of q , the parameter e is the local cone 
radius, i.e,, 

e s r s S tanq , 


and the equivalent length becomes 




tan^q 


S 

tan^ q dS = y . 
o 
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The method used to calculate the gas transport properties is based on 
the Lennard-Jones potential intermolecular force model as discussed in 
Ref. A-4. 


# Turbulent Convective Heat Transfer 

A similar method of analysis is used in this case as in the laminar case 
previously discussed and is presented in Ref. A*1 and A-5. The basic heat trans 
fer equation to be solved is Eq.(^-T)«, The Stanton number defined by Eq, (A-4), 
for turbulent flow becomes (Refs. A-1, A-2 and A-5). 


St = 0.0296 Re 

X 



(A-ijy 


Substitution of Eq.(A-ll)into Eq* (A-4), yields the enthalpy conduction as 


„ - 0»0296 
*h • p,2/3 


(PV^) 


0.8. ji 

' X 


0.2 
) . 


(A- 12) 


Substitution of Eq. (A- 12) into Eq. (A-1 ) and using the characteristic length 
(Ref. A-5) to account for nonuniform inviscid flow the turbulent flow convective 
heat transfer equation becomes, 


0.0296 .0.8 , „ , V 

q * (P <3^) - V • 

Pr ^ 


(A- 13) 


The starred parameters density, viscosity and Prandtl number are again 
calculated at a reference temperature corresponding to the Eckert reference 
enthalpy Eq. (A-5). The Eckert reference enthalpy is calculated as in the 
laminar case except that the recovery factor, r, is 


r = Pr 


1/3 


The recovery enthalpy is calculated using Eq. (A-2). 
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The stretched characteristic length, Xj is defined (Ref. A- 5^ by 

S 

*T = ‘-T 2 5 / »* V e dS. < A- 14) 

where S is the distance along the streamline. 

For a flat plate, the parameter e is set equal to a constant and thereby 
divided out, and for constant inviscid flow properties, the characteristic length 
X,p becomes equal to the flat plate length. 

For a cone with a cone half angle of q, the parameter e is equal to the 
cone local radius and is defined by 

€ = r = S tanq . 

For constant inviscid flow properties, the characteristic length X,p 
becomes 


X 


T " 


S'*25tanl*2% 


S 

f tan^’^^q dS = (|) S. 


o 
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